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A coordinate space approach, based on that used by Efimov, is applied to three-body systems 
with contact interactions between pairs of particles. In systems with nonzero orbital angular mo- 
mentum or with asymmetric spatial wave functions, the hyperradial equation contains a repulsive 
1/r 2 potential. The resulting wave functions are used in a renormalisation group analysis. This 
confirms Griesshammer's power counting for short-range three-body forces in these systems. The 
only exceptions are ones like the 4 S channel for three nucleons, where any derivatives needed in the 
interaction are found to be already counted by the scaling with the cut-off. 
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Effective field theories (EFT's) are now being widely applied to few-nucleon systems, see Refs. PJQ. The starting 
point is usually an organisation of the terms in an effective potential according to naive dimensional analysis (NDA) , 
as originally suggested by Weinberg Q . This classifies terms according to number of powers of low-energy scales they 
Oh' contain. In some cases, most notably iS-wave nucleon-nucleon scattering, the leading-order (LO) terms turn out to 
be unnaturally large. This is a consequence of low-energy bound or virtual states. It means that the LO terms need 
to be iterated to all orders in solving the Schrodinger or Lippmann-Schwinger equation |J, |5|, |6j . 

However it is now clear that NDA is not valid in all systems. There can be nonperturbative effects associated with 
^ i strong long-range potentials that significantly change the power counting for short-range interactions. This was first 
noted in the context of attractive three-body systems (such as three bosons, or the 2 S channel for three nucleons) 
where the leading three-body forces must be promoted to LO. More recently, failures of NDA have been observed in 
' repulsive three-body systems || , and for nucleon-nucleon scattering in spin-triplet waves 0, ^| . In the first example, 
£f) , short-range three-body forces in most low partial waves arc demoted to higher orders than naively expected; in the 
second, short-range forces are promoted to lower orders in P- and D- waves. In both cases, the terms in potentials 
O^N scale with noninteger anomalous dimensions, and so the standard classification of terms as LO, next-to- leading order 
(NLO) etc. is no longer convenient. 

In this note I examine repulsive three-body systems using the renormalisation group (RG) approach developed in 
Refs. |l2l ll'l Il4| . This provides an independent confirmation of results of Griesshammer for the power counting 
' in these systems 0. That work solved the Skorniakov-Ter-Martirosian equation ^jjj in momentum space, whereas 
here I work in coordinate-space following the approach developed by Efimov for attractive three-body systems [T^| 
' q ■ and recently extended by Gasaneo and Macek to cases with nonzero angular momentum I a l so comment on 

differences between the counting for derivative interactions in systems with strong long-range forces compared with 
the pure short-range case. This corrects the counting in Ref. [!j for the leading three-body force in the 4 S channel 
£^ | for three nucleons. 

•rH . If particles interact only through zero-range forces, then their wave functions satisfy the free Schrodinger equation, 
except where two of them coincide. The two-body forces can then be represented by boundary conditions at these 
points. These boundary conditions form the basis for Efimov's approach [l6| as well more recent work in Refs. [TtI 
liH Il9| . In particular Gasaneo and Macek have used this method to find solutions for systems with symmetric spatial 
wave functions. Here I generalise their results to cover asymmetric cases, such as the spin-quartet channels for three 
nucleons. It is convenient to work in hyperspherical coordinates since the boundary conditions are separable in the 
limit of infinite two-body scattering length. The resulting hyperradial equation then has the form of a free radial 
Schrodinger equation with a centrifugal- like 1/r 2 term whose strength is given by the hyperangular eigenvalue. This 
potential determines the form of three-body wave functions at small hyperradii, and hence it controls RG flow of the 
short-range three-body forces 0, Q . 

Sets of relative coordinates for three particles with equal masses are 



Xi=r fe -rj, y, ; = r t - - (r, +r fc ) , (I) 



where i,j, k are a cyclic permutation of 1,2,3, and I have used the traditional "odd-man-out" notation to label the 
sets. Hyperspherical coordinates can then by defined in terms of these as 




yf, on = arctan — . (2) 
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The three-body wave function can be decomposed into Faddeev components |20( as 
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*(r,n) = X] 



■i r 2 sin(2ai) 



fa(r, ai,Xi,yi), 



(3) 



where a factor of 1/ (xiUi) has been taken out to simplify the radial parts of the Hamiltonian. Away from the configu- 
rations where two particles coincide, each of these components satisfies a free Schrodinger equation. In hyperspherical 
coordinates, this has the form 
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Here Lj£, denotes the relative angular momentum of the pair jk, and Lj the angular momentum of particle i relative 
to that pair. If the problem is separable, we can write fa in the form 



0i(r,at,Xf,y<) = F i (r)ui(a i )Y l '. m >.(%)Y limi (y i ), 
where Fi(r) and Ui(aii) satisfy the ordinary differential equations 
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This hyperradial equation looks just like a free radial Schrodinger equation in two dimensions, with the hyperangular 
eigenvalue v 2 determining the strength of the centrifugal- like 1/r 2 term. In the cases of interest, where pairs of 
particles interact only in ^-waves, we can simplify the hyperangular equations by setting l[ = and l\ = I2 = h = I- 
The component fa is then independent of Xj . 

From the definition of the reduced Faddeev components fa in Eq. (J3J), they must vanish at Xi = 0. In hyperspherical 
coordinates these boundary conditions are 



«( r *I'*) 



0. 



(7) 



In the limit of infinite two-body scattering length, the logarithmic derivative of the reduced wave function must vanish 
whenever two particles coincide and so can interact via the two-body force [Tij Il7l ITsL Il9 | : 
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The points where x% = correspond to Xj = yj, y 3 = — iy.j and = — y$, yu = —hyi m the other relative coordinate 
systems. In terms of the hyperangular coordinates these are ay = = ^. The resulting boundary conditions on the 
Faddeev components are 



da* 



a ,=o sin(2aj) 
These are separable and lead to the hyperangular conditions 
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a fc =7r/3,y fc = -yi 



= 0. 



(9) 
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The symmetries of the spatial wave function can be used to simplify these conditions further. There are two cases 
of physical interest. The first is a spatial wave function that is symmetric under exchange of any pair of particles. In 
this case the Faddeev components are equal: 

fa = fa = fa = fa (11) 

This describes three identical bosons, or three fermions with different quantum numbers (spin, isospin, etc.) whose 
intrinsic state is completely antisymmetric. Most importantly for nuclear physics this corresponds to three nucleons 
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with total spin |. The second is where the spatial wave function is antisymmetric under exchange of one pair of 
particle, but symmetric under exchange of either of the others. The components are then related by 

01 =0, 2 = -0 3 = 0- (12) 

This describes three nucleons with total spin | . Generically we may write the boundary conditions in the form 



du 
do 



a=0 



+ A(-1)'^«(!)=0, (13) 



where A = +1 for completely symmetric spatial wave functions and A = — \ for cases with one antisymmetric pair. 

As noted by Gasaneo and Macek the hyperangular equation for I' = can be solved in terms of hypergeometric 
functions. After defining the new variable z = cos 2 a and writing u(z) = z( l+1 ^ 2 v(v), the equation takes the form of 
the hypergeometric equation |2l| . 



dv r , , ,d 



v 



with 



z(l - z)—r + [c - (a + b + l)z]— - abv = 0, (14) 
dz z dz 



l + l-v , l + l + v ,3 . ir , 

a= g ' b= 2 ' + 2' ( } 



The hyperangular eigenfunctions are thus 1 

fl + 1 — v l + l + v 3 9 \ 
u(a) = (cosa) + 2 F X I , ,1 + -; cos a I . (16) 

The corresponding radial solutions are just Bessel functions of order v and so we get 

0(r,a,y t ) = NJ v (pr) (cosa)' +1 2 F 1 ( ^ + 1 ~ V , l + 1 ^ t ' V ,l + |;cos 2 o^ Y lm {%). (17) 

These vanish at a — |, as required by the first boundary condition. The other condition, arising from the contact 
interactions at x± = 0, then provides an equation for the eigenvalues v 2 . The hypergeometric functions have the 
properties plj. 

— 2 F 1 (a 1 b,c;z) = ^- 2 F x {a + 1, b + 1, c+ 1; z), (18) 

and for z equal to or close to one 

, , . T(c)T(c - a-b) 

2 F 1 (a,b,c;l) = 1 1 r( ifc-a-6>0, 

1 (c — ajl (c — o) 

2 F 1 (a,6,c;z) ~ ^M£^+lz c ) if c _ a - & < 0, (19) 
Using these, Eq. (|13|l can be expressed in the form 

which matches Eq. (2.18) of Ref. 0, with the substitution of s by za It also agrees with the results of Ref. if A = 1. 
For symmetric systems (A = 1) with / = this is the equation first derived by Danilov [22], which has an imaginary 
solution for v. The corresponding hyperradial wave functions show oscillatory behaviour at small distances, and this 
is responsible for Thomas [2j| and Efimov 0] effects. For I > 1 or systems with one antisymmetric pair, the roots 



1 Note that I have to chosen t o wr ite these in the form that will make most direct contact with the results of Ref. rather than the 
equivalent form given in Ref. Il7l . Also, the factor of 2 in Eq. (13) of Ref. is incorrect and should be omitted. 
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of the equation are real and can be found in Table 2 of Ref . [!j . In these cases the 1 /r 2 potential is repulsive and so 
there is no Efimov effect. 

Having constructed the wave functions for the long-range forces in these systems, we can now use the methods of 
Ref. |13| to find the RG eigenvalues, which give the power counting for terms in the short-range potential. In fact for 
real values of v, we can immediately use the results in Eq. (54) of that paper if we multiply the hyperradial solutions 
by y/ir/ {2pr) to get functions that satisfy a three-dimensional radial Schrodinger equation. A term in the rescaled 
potential proportional to p 2n (n powers of the energy) varies with the cut-off A as A 2 ^ n+U ^ p 2n and so its RG eigenvalue 
is 

p = 2(n + v). (21) 

If we assign A-independent terms to LO in our expansion of the EFT, then p also labels the order of a term. The 
leading term in each channel is thus of order 2v, in agreement with the results in Table 3 of Ref. , except for the 
4 5* and Wigner-antisymmctric 2 S channels, where Griesshammer adds two extra powers of low-energy scales. 

The motivation for adding these two powers is the antisymmetry of the wave function in these channels which 
prevents all three particles from coinciding. As a result a pure <5-function interaction has no effect on them; at one 
with at least two derivatives would be needed. For any finite cut-off, however, a contact interaction becomes nonlocal 
and so can contribute. In Ref. 0, this happens implicitly through the momentum cut-off. In contrast, Ref. |l8j | does 
it explicitly by using a J-shell form for the short-distance interactions. This was done to ensure that the interaction 
has an effect even though the wave functions vanish as r — > as a result of the 1/r 2 potential. For example, the same 
RG analysis can be applied to two-body scattering with nonzero angular momentum L by setting v = L + i |l3j | . It 
shows that the leading short-distance interaction in this partial wave is of order 2L, as expected from the fact that 2L 
derivatives of a ^-function are needed to form a contact interaction that acts in this wave. Note that these derivatives 
are already counted by the RG eigenvalue, p = 2L + 1. 

The wave functions in two-body channels with nonzero angular momentum, or in three-body channels with v > — ^, 
satisfy a radial equation of the form 

Id 2 v 2 - i 

- T2 NO = -pfy + 1>. (22) 

r ar r z 

Acting on one of the wave functions where the long-range 1/r 2 interaction has been iterated to all orders, an interaction 
with two derivatives thus gives rise to two contributions. One is just proportional to two powers of the low-energy 
scale p, and so is two orders higher in the power counting. However the other, proportional to (1/r 2 )?/;, is of the same 
order as the term without those derivatives, since at small distances 1/r is not a low-energy scale. This second piece 
is absent if L — (or equivalently v = i). In that more familiar case, additional derivatives do indeed increase the 
order of the interaction. 

The bottom line is that any derivatives needed to construct appropriate interactions for the repulsive 1/r 2 potentials 
are already counted in the RG eigenvalue (or by the superficial degree of divergence in Ref. Q), without any need to 
add additional powers. But, apart from this rather minor amendment, the present analysis confirms Griesshammer's 
results for the power counting in repulsive three-body systems Ref. Q. The leading term in each channel has 
RG eigenvalue 2v, where v 2 in the hyperangular eigenvalue. These eigenvalues are not integers and so the usual 
classification of terms as NLO etc. in the EFT becomes inconvenient. In most cases this counting demotes short- 
distance three-body interactions to higher orders than predicted by NDA, although in some channels there is a small 
degree of promotion. 
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